Quantum Speed Limit With Forbidden Speed Intervals 
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Quantum mechanics imposes fundamental constraints known as quantum speed limits (QSLs) on 
the information processing speed of all quantum systems. Every QSL known to date comes from the 
restriction imposed on the evolution time between two quantum states through the value of a single 
system observable such as the mean energy relative to its ground state. So far these restrictions 
only place upper bounds on the information processing speed of a quantum system. Here I report 
QSLs each with permissible information processing speeds separated by forbidden speed intervals. 
They are found by a systematic and efficient procedure that takes the values of several compatible 
system observables into account simultaneously. This procedure generalizes almost all existing QSL 
proofs; and the new QSLs show a novel first order phase transition in the minimum evolution time. 

PACS numbers: 03.65.-w, 03.67.-a 



Introduction — Time is a valuable and often an ir- 
replaceable resource. Minimizing runtime is one of the 
most important driving forces behind hardware, software 
and computational complexity researches. Quantum me- 
chanics, as a fundamental law of nature, gives ultimate 
constraints known as QSLs on the runtime and hence in- 
formation processing speed of a computer, classical and 
quantum alike [l|-|3- So far, all QSLs are found by study- 
ing the evolution of quantum systems given the value of a 
certain observable. Interestingly, they all give evolution 
time constraints in the form r > g^(F, v), where r is the 
time needed to evolve from one normalized state to 
another normalized state \<p) with fidelity F = | (ip | ip) | 2 , 
and v is the value of a certain observable O of the sys- 
tem The significance of these results is that know- 
ing partial information of the system as little as the value 
of just one observable of the system is enough to deduce 
its maximum possible quantum information processing 
speed. For example, using the energy standard deviation 
AE as the observable, the corresponding QSL, namely, 
t > Kcos~ l (V~F)/AE, is the famous time-energy uncer- 
tainty relation [5rl9j. 

The more information about the quantum system is 
given, the more stringent the QSL one should be able 
to obtain. So it is instructive to investigate what kind 
of QSLs one can deduce when more than one compat- 
ible observable of the system are given. However, very 
limited progress has been made along this line. All rele- 
vant works to date simply consider the constraints that 
logically follow from the QSLs for each of the observ- 
ables [1, 14 1 instead of analyzing the restrictions due to 
these compatible observables holistically. 

Here I report a powerful method to st udy QSLs by 
extending earlier proof techniques [U, EH, EH- Besides 
finding QSLs when several compatible observables of the 
system are given, this method also provides a unified way 
to show all known QSLs. Through these new QSLs, I find 
that the minimum possible evolution time can exhibit a 
new first-order phase transition with fidelity F being the 



order parameter. Finally, I will define and briefly discuss 
the reverse problem of QSL construction. 

Construction of an auxiliary inequality — I write 
the time-independent Hamiltonian H in the diagonal 
form Y). Ej\Ej)(Ej\. Surely, any normalized quantum 
state \ip) can be expressed in the form Y) , di\Ef) with 
J2j \a 3 \ 2 = 1. Suppose \tp) — e~ lHT / h \ip) is the state ob- 
tained by evolving \ip) by H for a time r > 0. Then, the 
fidelity F between and \<p) obeys 
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for any real-valued 6. Here I have used the fact that 
the magnitude of a complex number is greater than or 
equal to its real part to arrive at the above inequality. 
Actually, inequality ([TJ is an extension of those used in 
Refs. 

Let p(x) be a function satisfying p(x) < cos(x — 9) 
whenever x > 0. Then, 
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provided that Ej > for all j. The following section 
shows that this p{x) can be chosen to be a polynomial-like 
function in the form Y^k=o c k xSk with s > efficiently. 
(This is a mo re g eneral choice for p(x) than in all previous 
studies (l-EEflHEi-) 

Proof of the existence of efficiently computable p(x) 
- It suffices to show the existence of a polynomial 
q(x) < cos(x 1 / s — ff) = f(x) whenever x > 0. In 
fact, q(x) exists even if I further demand that it meets 
f(x) at finitely many distinct non-negative points, say, 
xts. I do this by considering the Hermite interpolat- 
ing polynomial q(x) that satisfies cf k \xi) = f^ k \x) for 
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k = 0, 1, ... , 2j( — 1 for some h € Z + . This polynomial 
can be constructed efficiently 15[. Surely, f(x) — q{x) = 
ag(x — xi) 2je + 0((x — xg) 2 ^ e+1 ) locally around each X£. 
For randomly chosen xg's and j^s, all the ags are non- 
zero almost surely. And in the singular case in which 
some of the ag equals 0, I simply randomly choose an ex- 
tra distinct point x a and demand further that q(x) obeys 
q {k) {x a ) = f (k) {x a ) for k = 0, l,...,2j a - 1 for some 
randomly chosen positive integer j a . Then, the modified 
q(x) locally agrees with f(x) up to an even power of x—Xf, 
at each of the interpolating point xi almost surely. 

Recall that there are efficient and stable algorithms 
to find all the real roots of a polynomial [l6j]. Apply 
one such algorithm to find the largest real root x u of 
the equation q~(x) = — 1, I can bound the non- negative 
roots of q(x) = f(x) to the interval [0, x u ]. Since f(x) = 
cos^x 1 ^ — 9) is a smooth function of bounded variation in 
[0, x u ], I may use interval arithmetic to efficiently find all 
the sub-intervals of [0, x u ], if any, in which f(x) — q(x) is 
negative [T3| • Actually, there are at most a finite number 
of these sub-intervals; and they are present if and only if 

• /(*) - q(x) = a e (x - x e )*** + 0((s - x t f h+l ) in 
the neighborhood of xe with ag < 0; or 

• for two consecutive distinct roots x\ and xi of 
fix) — q(x), there are x± < yi < y% < X2 such 
that f(x) — q~(x) < for all x € (2/1,2/2)- 

In the first case, I may bring f{x) — q(x) up above zero 
by adding a term bg(x — xg) 2 ^ 1 with bg > at,. Whereas 
in the second case, this can be done by adding a term 
in the form V U e (x - x t ) 2(j(+K ^ with K t = 1 if a t > 
and Kg — if < 0. Note that this additional term 
does not affect the local behavior of f(x) — q(x) around 
those xis with ag > 0. Since there are only a finite 
number of such sub-intervals, I can efficiently find b > 
such that f(x) — q(x) > for all x > where q(x) = 

q(x) - bUt(x ~ xe) 2Ut+Kt) ■ 

Last but not least, I remark that since what one re- 
ally need is q(x) < f(x) for x > 0. So, whenever 
Xf, = 0, namely, the boundary point, there is no need 
to demand f(x) — q(x) — a t x 2 ^ 1 + 0(x 2je+1 ) for x suffi- 
ciently close to 0. Suppose f(x) — q(x) = ai J + 0(a; : ' +1 ) 
in the neighborhood of x — for some positive inte- 
ger j and a 7^ 0. Then, the transformation q(x) — 
q(x) - bx j+K Y\' e (x - xg) 2( - je+K ^ for a sufficiently large 
b > will do. Here n = 1 if a > and k = otherwise. 
Besides, the primed product is over all xg 7^ 0. 

Construction of the QSL from p{x) — Substituting this 
polynomial-like p(x) into inequality ([2]), I conclude that 
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the energy of the system. Furthermore, in the case 
of 9 = 0, I may rewrite inequality ([T]) as \AF > 
J2j \ a j\ 2 cos(\Ej\T/H). So the above arguments lead to 
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irrespective of the signs of E/s, where = 
J2j \ a j\ 2 \Ej\ r - I remark that inequality ^ becomes an 
equality if and only if e~ l6 (ip\(p) is real and non-negative 
together with cos(Ejr/h — 9) = p(EjT/H) for all j with 
ocj ^ 0. 

Consequently, suppose the values of the compatible 
(time-independent) observables of the system (E k ) (or 
(\E\ k )) are known for k = 1,2, ...,n. Then, given a 
fixed fidelity F £ [0,1], the required evolution time must 
satisfy inequality @. Since there are efficient numer- 
ical algorithms to find real roots of a polynomial equa- 
tion [18[ , I can find the permissible intervals for r readily. 
As the reference energy level has no physical meaning, I 
may tighten the permissible region for r by taking the in- 
tersection over all the permissible intervals obtained by 
replacing (E k ) by ((E + a) k ) (or (|£| fe ) by (\E + a\ k )) for 
all a and k — a trick first used in Ref . [12| . Finally, I may 
further strengthen the bound by taking the intersection 
over the permissible regions for r obtained by all degree 
< n polynomials p(x) < cos(x — 9) for x > 0. There 
is no known efficient way to perform this very last task, 
however. 

Recovering all existing QSLs — The above procedure, 
even without the final step, is already powerful enough to 
prove all the known QSLs. I choose p(x) to be the func- 
tion 1 — ax b which meets the curve cos(a; — 9) at two dis- 
tinct points for x 7^ 0, namely, at x = and x c such that 
p(x) actually touches the curve cos(a; — 9) tangentially 
at the latter point. Then, in the event that 9^0, in- 
equality (|3"a|) gives the Margolus-Levitin bound [2|-|4[ and 
its generalization [ll|; whereas in the event that 9 = 0, 
inequality (13bl) becomes the Chau bound [13] and its gen- 
eralization [13| . 

To recover the time-energy uncertainty relation, I 
use the inequality cos a; > 1 — x 2 /2. From inequal- 
ity (3l|, I get the bound </F > 1 - {\E\ 2 ) [t / h) 2 / 2 = 
1 — (E 2 )(t/H) 2 /2. This bound can be optimized by choos- 
ing the reference energy level to be the average energy of 
the system. The result is V~F > 1 - (tAE/H) 2 /2 pro- 
vided that the system evolves under a time-independent 
Hamiltonian. Now I consider evolving the system for an 
infinitesimal time At. The constraint set by the above 
inequality becomes VF > cos(ArA J B/ft) + 0((Ar) 4 ) = 
cos(AtAE/H + 0((At) 3 )). Hence, the corresponding 
infinitesimal change in Bures angle must obey Ai? = 



whenever Ej > for all j, where (E r ) = J2j \ a j\ 
denotes the expectation value of the rth moment of 



cos _1 (VF) < AtAE/H + 0((At) 3 ). Since Bures an- 
gle is a metric (l9j . by integrating over a finite time, 
I conclude that for a time-dependent Hamiltonian, the 
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(a) (b) 




FIG. 1. (a) The time evolution of root fidelity \J~F for [<p e (0)) under the Hamiltonian H e . (b) The polynomial p e {x) is a very 
good approximation of cosx for < x < 22. (c) The curves cos(a;) and p^ 2 \x) intersect at four distinct points only, (d) A 
plot of the R.H.S. of inequality (3b]) with p (») = p e {x) and (\E\ 2k ) = 4h 2k [1 + (ll/5) 2fc ]/15 for k = 1, 2, . . . , 17. 



evolution time r obeys JT AE 1 dr > hd — hcos~ x (V~F), 
which is the time-energy uncertainty relation for time- 
dependent Hamiltonian. If the Hamiltonian is time- 
independent, the above expression becomes the famous 
inequality r > h cos^ 1 (V r F) / 'AE. 

Existence of QSLs with forbidden speed intervals — By 
picking the function p(x) with n > 1, I have the surpris- 
ing situation that the permissible evolution time r given 
by inequality ([3]) is in general separated by forbidden 
time intervals. This is not completely unexpected be- 
cause unlike all previous QSLs, here the quantum system 
is constrained by more than one compatible observables. 

One may question the genuineness of these forbidden 
time intervals as some of the apparently permissible time 
intervals are illusory because they could be the result 
of poorly chosen 9 and p{x). Nevertheless, the example 
below shows the contrary. 

Consider the initial state |^ e (0)) = [V7\0) + V2(\-h) + 
\K) + l-im/5) + |HV5))]/Vl5 evolving under the 
time-independent Hamiltonian H e = J2j=o ±i ±11/5 
HEj\hEj)(hEj\. FigOJi depicts the time evolution curve 
for the root fidelity >/F = \(tp e (Q)\ip e (T))\, show- 
ing that the first time for \/~F to reach \JF C \ = 
and \/F C 2 ~ 0.0682 are at r = t c i ~ 9.693 and 
t C 2 « 4.110, respectively. I choose p(x) = p e (x) 



to be the Hermite interpolating polynomial satisfying 
the following constraints: p^\o) = d 3 p e (0) / dx^ = 
cos^(O) for j = 0,1,2, p ( e 3 \±T cl ) = cos( j ')(±t c i) for 
j = 0,1,2,3, p { e'\±T c2 ) = cos«(±r c2 ) for j = 0,1, 
pi j) (±llr cl /5) = cos«(±llr c i/5) for j = 0,1,2,3,4,5, 
and ^(±11^/5) = cos( J )(±llr c2 /5) for j = 0,1,2,3. 
By construction, p e (%) — cos a; at x € S e = 
{0, ±t c i, ±t c2 , ±llr c i/5, ±llr c2 /5}. By construction, 
p e (x) is even and of degree 34. Fig. [lb shows that p e (x) is 
a very good approximation to cos a; for < x < 22. More 
importantly, I show in the next section that p e (x) < cosx 
for all real x. 

Proof ofp e (x) < cosx — It is straightforward to check 
that p e (x) > cosx in the neighborhood of x g S e . To 
show that this is also true for all x £ K, I need the 
following lemma. 

Lemma 1. Let f{x): R — > R be a real-valued differen- 
tiable function with exactly n real roots counted by multi- 
plicity. Then, f'{x) has at least n — 1 real roots counted 
by multiplicity. 

Proof. The lemma is a simple consequence of the follow- 
ing two facts. First, if x\, x 2 are two distinct consecutive 
roots of /, then Rolle's theorem implies that there is a 
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root x € [xi,X2] for /'. Second, suppose x\ is a multiple 
root of / of multiplicity k > 1, then clearly x% is a root 
of /' of multiplicity k — 1. □ 

By construction, x = is a root of multiplicity 3 for the 
even function g(x) — p e (x) — cosx. Similarly, by counting 
the multiplicity of roots of g(x) at a; € S e \{0}, I conclude 
that g(x) has at least 35 real roots. Suppose it had more 
than 35 such roots, then the even function g(x) should 
have at least two more real roots — one positive, one 
negative. By Lemma [1] g^ 2 \x) — p^ 2 '(x) — coax would 
have at least 35 + 2 — 32 = 5 real roots. However, a plot of 



cos x and the quadratic function pi ' (x) in Fig. [TJ; shows 
that g( 32 '(x) only has 4 roots counting by multiplicity in 
the range x G [—8,8]; and g^ 32 \x) does not have any 
real root outside this range as pi (x) < —1 for |x| > 8. 
Thus, S e is the set of all real roots of g(x). Since g(x) < 
in the neighborhood of these roots, the continuity of g 
implies that g(x) < for all x € M. 

Actually, this method can be adapted to show p(x) < 
cos(x — 9) for all x > for a variety of polynomial p(x) 
constructed out of Hermite interpolation. 

QSL and a first- order phase transition - Since 
Pe(x) < cosx for all x > 0, it induces the following QSL: 



VF< 1 - 5.0000 x 10 1 < | ^1 2 > + 4.1667 x 1Q- 2 (\E\ 4 ) - 1.3889 x 10~ 3 ( 



+ 2.4802 x 10- 5 (|£| 8 ) Q) 8 - 2.7557 x 10- 7 (|£| 10 ) (^j™ + 2.0876 x 10~ 9 (|£;| 12 ) Q 



I ' y 1.5585 x lO-^d^l 18 ) ( T - 



1.1469 x 10- n (|£| 14 ) (~) + 4.7766 x 10" 14 (|^ 
4.0798 x 10 19 O^l 20 ) (-J - 8.6954 x 10" 22 (|S 

2.1160 x 10- 27 (|S| 26 ) (^) 26 + 2.2929 x 1 30 < 2S > {^j* - 1.79 55 x 1 33 < | 1 30 > (^) 



22 



1.5 1 25 x 10 - 24 (|£f 4 ) 



18 



24 



:S(> 



19531 x 10~ (LE 



^ I T V 2 2.1140 x 10- 40 (|£l 34 \ n " !i 



(4) 



In particular, by putting {\E\ 2k ) = (E 2k ) = Ah 2k [l + 
(ll/5) 2fe ]/15 for k = 1, 2, . . . , 17, this construction leads 
to a QSL which gives a tight lower bound for the evolu- 
tion time in the cases of VF = V-^ci and VFc2- Fig. QJi 
depicts that for root fidelity v^F = y/Fci — 0, the corre- 
sponding QSL is r e [t c1 , 10.138] U [10.248, +oo). Com- 
bined with the example of the evolution of |^ e (0)), I 
conclude that whenever a quantum state with (|£'| 2fc ) 
is equal to the value given in the previous paragraph 
for k — 1,2,. ..,17, the minimum evolution time r m j n 
for it to evolve to another state of root fidelity is t c \. 
By gradually increasing yF, the allowable region for 
r increases and r m i n < t c \. Most importantly, by in- 
creasing \f~F to \JF C 2 « 0.0682, the allowable r becomes 
{t"c2} U [9.519, +oo) with r m ; n = r C 2. That is, a new per- 
missible time interval appears and a genuine forbidden 
evolution interval (t C 2, 9.519) is formed. Besides, r m - m 
shows first-order phase transition at y/F = \TFa- This 
is the first QSL that captures this type of phase transi- 
tion. Further significance of this result is reported in the 
Supplementary Material [2pj ] . 

The reverse problem — Note that the above method to 
construct a QSL with genuine forbidden speed intervals 
is generic. In fact, it brings us to the following reverse 
problem, which has never been studied before. Given an 



initial state, an Hamiltonian and a required fidelity F, 
is it possible to find a QSL whose minimum permissible 
evolution time equals the actual evolution time needed? 
By modifying the proof of the existence of p(x), I show 
in the Supplementary Material that the answer is affir- 
mative in finite-dimensional Hilbert space [2(ij |. 

Outlook — In addition to showing that QSL can have 
forbidden speed intervals, this work opens up a more 
general research direction on the tradeoff between the 
amount of partial information describing a quantum sys- 
tem and the constraints on its information processing 
capability in which a lot of works can be done. 

I thank F. K. Chow, C.-H. F. Fung and C. Y. Wong for 
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by the RGC Grants HKU 700712P and HKU8/CRF/11G 
of the Hong Kong SAR Government. 
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Appendix B: Existence Of A QSL For The Reverse 
Problem For Finite-Dimensional Hilbert Space 
Systems 

Denote the state at time r under the evolution of 
the time-independent Hamiltonian H in a d-dimensional 
Hilbert space by \(f(r)) with \(p(0)) = J2j=i a j\-^j) ■ 
Surely, the root fidelity at time r is given by the con- 
tinuous function 



Appendix A: Further Significance Of The QSL 
Associated With The Degree 34 Polynomial p e (x) 
Reported In The Main Text 



Recall that for fidelity F — 0, the QSL reported in 
inequality (4) in the main text leads to a minimum evo- 
lution time of T c i for (\E\ 2k ) = Ah 2k [l + (ll/5) 2fc ]/15 
for k — 1,2,..., 17. Further, this bound is tight for 
it can be achieved by the state |<^ e (0)) under the evo- 
lution of the Hamiltonian H e . Note that the highest 
and lowest energy eigenvector components of |<yj e (0)) are 
| ± HH/b). Hence, the phase angle difference \ rotated 
during the time r c i between these two components equals 
22t c i/5 > 2tt. That is, the relative phase angle between 
two components has to rotate more than one complete 
circle in order to evolve \ip e (0)) to its orthogonal com- 
plement. This is a new situation for these relative phase 
angles rotated in all known QSLs to date [lMff are at most 
2tt. The implication is that for states obeying the above 
constraints on its various moments of energy, they can- 
not evolve to their orthogonal complement without some 
time of "time wastage" as some of the relative phase an- 
gle change must be greater than a complete circle. 

I also remark that this is the first tight QSL with the 
property that the "magic state" saturating this QSL in 
the case of F — has to be at least four-dimensional. 
The corresponding "magic states" for all previous QSLs 
are at most three-dimensional [IrS]. The reason why 
the "magic state" saturating this QSL is at least four- 
dimensional is as follows. From the discussion on the 
conditions for equality of inequality (3b) in the main text 
and the construction of p e that leads to the QSL, the 
magic state \ip(0)), expressed in the energy eigenbasis, 
must be in the form a |0)+Q!i| J B>+Q!2|-^)+a3|ll J B/5) + 
— ll-E/5) with E > 0. Furthermore, the evolution 
time to an orthogonal state equals T c \hjE. For the given 
constraints in (|i?| 2fc )'s, I arrive at |a | 2 = 7/15, \ai\ 2 + 



VF(r) = |M0)|^(r))H^K| 2 cos 



where 9(t) is the argument of the complex number 
(t^(0) |^(7-)) . Although 9 can only be determined mod- 
ulo 2tt, I may uniquely fix it by the integral curve describ- 
ing the time evolution of the argument of ((p(0)\(p(T)) 
with initial condition 6(0) = provided that r is less 
than or equal to the first time when y/F = 0. And from 
now on, I assume 9(t) to be this smooth integral curve. 

I first write down several properties of the function 
y/~F(r). Denote the first time when y/F(r) reaches a 
certain fixed value \/Fq 6 [0, 1] by T min . Suppose fur- 
ther that r m i n is finite. I define rt urn and eo as follows. 
Suppose y/F{r) is a decreasing function in [0, r m i n ], then 
I set T tum = and eo = 1. Otherwise, since y/F is 
continuous and differentiable provided that y/F > 0, 
there is a turning point in [0,r m i n ). Denote the turn- 
ing point in [0,r m i n ) closest to T m - ln by T turn . Then, I 
set e = mm Te[o. Ttu,-„] vF(t) — V^o- It is well-defined 
because the minimum exists owning to the continuity of 
y/F ; and it is positive for otherwise r m i n will not be the 
first time when the root fidelity reaches y/Fo. Note that 
no matter whether there is a turning point for y/F in 
[0, r m i n ] or not, y/F(r) is decreasing in [r turn , r m i n ]. More 
importantly, for < r < r m i n , y/F(r) < %/^o + eo only in 
[Ttum ) Ttnin]. That is to say, the function y/F(r) is one- 
one in the domain [r tU rn' , T min ] and range [y/Fo, y/Fo+eo], 
where r turn ' is the closest point to r m ; n in [0,r m i n ) with 
V^F(Tturn') = y/Fo + eo- Last but not least, for a suffi- 
ciently small 5 > 0, 



VF(T)-y/F~o = 0((T min -Tf) 



a 2 \ = |Q£ 3 | +KI =4/15. Consider the imaginary for r € (r min - 5, r min ) for some /3 > 0. 



part of (%p(0)\ip(T c ih/ E)) , I have (|a 2 | 2 — |ai| 2 ) sinr c i + 
( |o:4 1 2 — |a3 1 2 ) sin(llr c i/5) = 0. Hence, at most one of 
the ay's can be zero. Thus, the "magic state" is at least 
four-dimensional. 



Next, I consider the QSL induced by a polynomial p(x). 
Using the idea in the main text, suppose cos a; > p(x) = 
YZ=o Ck3;k for a11 x - x min = min,- mxa Te [ Q>Tm .^[EjT/H- 



2 



(r)]. Then, I have a QSL in the form of an inequality 
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d n 
n k 



F ^ >;k«il a p(^ I -«(r)) 



1 fe 



-EE 

fc=0 1=0 
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(B.l) 



whenever < r < r m j n . 

Now, I consider the set of polynomials 5n i7 with the 
properties that p(x) G S n ~ if and only if 

• degp(x) < n; 

• p{x) < cos x whenever x > 

• p(xj) = cosx,- for all j, where Xj = EjT m i n /h — 
0(T min ); and 

• cos a; — p(x) = 0((x — Xj) 1 ) in the neighborhood of 
Xj for all j. 

Using the Hermite interpolating polynomial construction 
in the main text, I know that for each fixed 7 > 1, the 
set S nn is non-empty provided that n is sufficiently large. 
Clearly, <S„ j7 is a convex set. In addition, it is easy to see 
that the functional 



G[/i,/ 2 ] 



\fi(x) - fr(x)\ 



is convex where x max = maxj [EjT m i n /h — 6*(r m i n )]. 

Note that for a fixed 7 > 1, there is a se- 
quence of polynomials p n (x) G iS n , 7 such that 
lim„_>. +00 G[p n (x), cos x] = 0. In fact, each p n (x) can 
be chosen to be the optimal degree < n polynomial in 
that minimizes the functional G[cosa;, •] via convex 
optimization Q. 

With the above background preparation, I am ready 
to prove the existence of a polynomial p{x) that solves 
the reverse problem. That is, the QSL induced by p(x) in 
inequality (jB.ll) has the property that the smallest non- 
negative time at which the R.H.S. of this inequality is 
\/Fq occurs when r = T m - m provided that (_E fc )'s are set 
to the kth moment of the energy of the state \ip(0)). 

I choose a sufficiently small < S < r min — T turn such 
that 

VF(t) - V^b > C(Tmin - Tf 

for all t G (T m ; n — S, T min ] , where ( > 0. For this 6 > 0, I 
can find sufficiently large 7 and n such that <S n , 7 is non- 
empty and the p(x) G 5„ i7 that minimizes the functional 
G[cosir, •] obeys G[cosx,p(x)] < e where 



S)- V^bG(0,e ]. 



I claim that the QSL induced by this p(x) solves 
the reverse problem. This is because by my construc- 
tion, for r = r min , the R.H.S. of inequality (|B.1[) equals 



E 



1 cosx. 



= y/F\ ) = Vf(i 



,„-,, 2 p(Xj) = J2j \0tj\ COSXj = VTo = V^^nin). 

I proceed to consider the case of r G [0, r m j n — 6]. 
cosx — p(x) < G[cos x,p(x)] < e for all x > x m \ 
conclude that 



As 
I 



v / F(r) = 



E 

j 



' cos 



Hence, 



Ei 



< e 



j 



(r)) > VF(r) 
= \fFh~. 



In other words, the R.H.S. of inequality (jB.ljl 
than y/Fo in the time interval [0, r m i n — S]. 

Finally, I consider the case of r G (r m j n — 6, r n 
have 

+ C(Tmin - rf 

=Ek^(^-*( t ))+°(( a ^)) 7 ) 



greater 



where Aj(r) = £y(r min .-T)/^+0(r)-( 
is smooth in this time interval, Aj(r) 
for some 7' > 0. Therefore, 



(r min ). Since 9{t) 
= 0((r min - t)t') 



E 



| 2 p( ^ - fl(r)) 



^ + C(r mi „ - rf + 0((r min - r)™'). 



Since £ > 0, by picking a sufficiently small 6 > and 
a sufficiently large 7 so that 77' > (3 (and a sufficient 
large n so that 5„ i7 is non-empty), the R.H.S. of inequal- 
ity (|B.1|) is greater than y/Fo in this time interval. 

To summarize, the QSL induced by any p(x) in this 
S nt /3 is a solution of the reverse problem. This completes 
the proof of my claim. 

Lastly, let me make the following remark. Suppose 
< T\ < r 2 < ... < Tg are t distinct numbers with 
V~F(Tj) — y/To for all j, where y/F(r) is the root fidelity 
between |v(r)) and |<£>(0)) under the action of a time- 
independent Hamiltonian H. Then, it is not difficult 
to adapt the above procedure to construct a polynomial 
whose induced QSL has the properties that 

• the induced QSL is an equality at times t\ , . . . , m 
provided that y[~F = y/Fo and (E k ) is the A;th mo- 
ment of the average energy of the state \ip(0)); 



3 



• the induced QSL is a strict inequality at time r € 
[0,Ti] provided that Vf(t) > F , 

The proof is left to the interested readers. 
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